AO-A167  S9C  DURABILITY  UNDER  REPEATED  BUCKLING  OF  STIFFENED  SHEAR 
PANELSCU)  TECHNION  -  ISRAEL  INST  OF  TECH  HAIFA  DEPT  OF 
AERONAUTICAL  ENGINEERING  H  KOLLET  ET  AL.  HAR  83 
UNCLASSIFIED  TAE-589  F/G  13/13 


AD-A167  596 


puRuni.iTY  paper  kepi- atfd  bpcflina;  of 
stiffs  \::d  mii  \k  panels 


by 


M.  kollet 
i.  Weller 
A.  I  ib:ii 
.  I.  S  i  urc  i 


Scientific  Report  Ho.  S 


r  ‘ **•  4!  ^  **  ■  -y  • .  ’  ,  >- 


...  • 

i  .c 

. *  ■■  H 

M  *«  t  f  sj  I  7  .  *  3*  1  . 


V 

'’'‘■iS* 

o, 


r 


Approved  for  public  release;  distribution  unlimited 


Prepared  for 


CD 


(  ,_> 
f. «. 


structures  and  Pynnraics  Pi  vision,  \ir  I erce  Wright  Aeronautical  Laboratories 

and 


Pur  >pean  Office  of  Aerospace  Re  search  and  development,  London  England 


March  1983 


DURABILITY  UNDER  REPEATED  BUCKLING  OF 
STIFFENED  SHEAR  PANELS 

by 

M,  Kollet 
T,  Weller 
A  Libai 
J,  Singer 


Department  of  Aeronautical  Engineering 
Technion  -  Israel  Institute  of  Technology 
Haifa,  Israel 


The  research  reported  in  this  document  has  been  sponsored  by  the  Air  Force 
Office  of  Scientific  Research,  United  States  Air  Force,  contract  AF0SR-81-0016 

Distribution  of  this  document  is  unlimited 


-  ii  - 

TABLE  OF  CONTENTS 

ABSTRACT  i 

TABLE  OF  CONTENTS  ii 

NOTATION  iii 

LIST  OF  TABLES  iv 

LIST  OF  FIGURES  v 

L.  INTRODUCTION  1 

2.  TEST  SETUP  5  PROCEDURES  8 

3.  TEST  SPECIMENS  12 

4.  NUMERICAL  STUDIES  14 

5.  TEST  RESULTS  -  DISCUSSION  §  ANALYSIS  18 

6.  CONCLUSIONS  AND  RECOMMANDATIONS  FOR  FUTURE 

WORK  23 

ACKNOWLEDGEMENTS  25 

REFERENCES  26 


TABLES  AND  FIGURES 


29 


ABSTRACT 


An  experimental  investigation  was  carried  out  on  the  durability 
of  thin  webbed  beams  subjected  to  repeated  shear  bucklings.  Although 
the  repeated  loading  exceeded  the  initial  buckling  load  several  times, 
it  was  in  the  range  of  25%  to  50%  of  the  ultimate  loads,  which  is  well 
within  the  flight  evelope  of  many  airframes.  Initiation  and  propaga¬ 
tion  of  fatigue  cracks  up  to  complete  failure  was  followed  during  the 
tests.  An  empirical  life  prediction  formula  based  on  this  test  series 
was  developed  and  will  be  updated  by  further  testing. 

The  experimental  investigation  was  supplemented  by  numerical 
analysis  using  the  STAGS  computer  program. 

..  t/ 

The  number  of  load  cycles  to  crack  initiation,  ,  and  the  fatigue 
failure  N^,  ranged  in  this  test  series  from  about  20,000  cycles  to  about 
200,000  cycles ,  .  ,  . 

It  appears  that  the  most  significant  parameter  affecting  Ni  was 
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the  ratio  of  ultimate  to  cycling  loads.  Tv -y  •T'  -  K ■  - 
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correlation  parameter 
amplitude  of  assumed  imperfection 
distance  between  uprights  (ss  Fig.  2) 
overall  height  of  Wagner  beam  (see  Fig.  2) 

area  moment  of  inertia  of  flange  and  upright  cross  sections  respectively 

torsional  moment  of  inertia  of  flange  and  upright  cross  sections 
respectively 

number  of  cycles  to  crack  initiation 
number  of  cycles  to  failure  in  fatigue 
web  thickness 
external  shear  force 

experimental  external  shear  force  at  buckling 
calculated  external  shear  force  at  buckling 
experimental  external  shear  force  at  ultimate  failure 
calculated  external  shear  force  at  ultimate  failure 
external  shear  force  at  which  local  yielding  first  occurs 
transverse  deflection  of  beam 
initial  imperfection 
correlation  parameters 
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1,  INTRODUCTION 

Stiffened  plates  and  panels  can  withstand  shear  loads  far  in  excess 
of  their  initial  buckling  loads  before  they  fail.  This  stable  postbuck- 
ling  behavior  has  been  an  important  factor  in  many  aeronautical  engineering 
designs,  especially  where  strength  is  the  primary  design  criteria  and  redu¬ 
ction  of  structural  weight  is  at  a  premium.  The  structural  model  which 
is  commonly  used  for  the  analysis  of  this  phenomenon  is  the  theory  of 
"diagonal  tension",  first  introduced  by  Wagner  [21  for  the  case  of  rigid 
stiffeners.  It  was  further  developed  by  Kuhn  and  his  associates  to  include 
the  effects  of  the  tension  field  on  the  compressive  load  carrying  capabi¬ 
lity  of  the  plate.  Also  included  were  the  effects  of  the  in-plane  flexibi¬ 
lity  of  the  stiffeners  and  the  participation  of  the  plate  in  the  loads  car¬ 
ried  by  the  stiffeners.  This  resulted  in  the  theory  of  "incomplete  diago¬ 
nal  tension"  which  is  summarized  in  [3-6] . 

During  the  fifties  and  sixties,  utilization  of  the  stable  postbuckling 
behavior  of  plates  in  shear  was  extended  to  civil  engineering  applications, 
mainly  for  plate  girder  design.  Some  results  of  these  studies  are  summari¬ 
zed  in  [1],  where  an  extensive  list  of  references  is  also  supplied.  In  most 
civil  engineering  designs  the  flanges  of  the  plate  girder  are  relatively 
flexible  and  the  mode  of  ultimate  failure  can  be  flange  collapse  rather  than 
web  tear.  Also,  the  depth  of  the  wave  pattern  in  the  web  is  affected  by  the 
flange  flexibility.  The  criterion  introduced  by  Rockey[7]  can  be  used  to 
estimate  the  importance  of  this  effect.  See  also  equation  (1)  of  [1]. 


The  designs  based  on  criteria  developed  by  tension  field  theory  and 
by  the  civil  engineering  criteria  of  flange  failure^ considered  the  static 
ultimate  load  carrying  capability  of  the  web-stiffener  system  in  shear  as 
the  prime  factor  to  be  taken  into  account.  Therefore,  all  the  analyses  were 
aimed  at  the  determination  of  the  ultimate  shear  load,  Details  of  the 

stress  distribution  at  the  working  loads  were  not  supplied  and  were  not  con¬ 
sidered  to  be  worth  the  effort  of  obtaining  them,  provided  that  the  ultimate 
margin  of  safety  was  adequate.  In  any  case,  these  details  were  difficult 
to  calculate  and  even  now,  with  the  extensive  capability  of  computational 
facilities,  "a  full-fledged  nonlinear  analysis  that  utilizes  the  necessary 
refined  grid  to  model  the  actual  structure  is  prohibitively  expensive  and 
therefore  unrealistic  from  a  computational  and  design  point  of  view [8]” 

Two  additional  factors  have  been  added  during  recent  years  which  have 
led  to  the  need  for  reappraisal  of  the  subject.  The  first  is  the  appearance 
of  new  materials  for  the  design  of  primary  parts  of  aircraft.  Among  these, 
higher  strength  metals  and  composites  are  the  more  noteworthy.  The  second 
(and  of  no  lesser  importance)  factor  is  the  steady  increase  in  the  number 
and  relative  levels  of  repeated  load  applications  during  the  expected  life¬ 
time  of  the  aircraft,  as  a  result  of  more  optimal  designs  and  of  better 
utilization  and  extended  life  of  airplane  fleets. 

The  main  concern  is  that  parts  which  are  designed  to  be  subjected  to 
many  repeated  bucklings  during  their  lifetimes  may  fail  prematurely  in  fa¬ 
tigue.  The  danger  increases  with  the  number  of  bucklings,  the  level  of  the 
shearing  loads  above  the  initial  buckling  load,  and  the  susceptibility  of  the 
material  to  local  effects  caused  by  the  high  localized  bending  stresses 


which  are  superimposed  on  the  in-plane  stresses.  The  problem  also  came  up 
in  civil  engineering  applications  [9-11]. 

This  concern,  supplemented  by  recent  service  failures  which  could  be 
attributed  to  repeated  bucklings,  provided  the  general  motivation  for  the 
present  study.  The  results  reported  here  are  confined  to  metal  beams,  and 
the  material  utilized  is  Aluminum  2024T3,  which  is  one  of  the  standard  ma¬ 
terials  in  aircraft  manufacturing.  Extension  to  composite  beams  will  be 
made  at  a  later  stage. 

Direct  literature  on  the  durability  and  fatigue  failures  of  stiffened 
shear  panels  subjected  to  repeated  bucklings,  is  not  extensive.  Reference 
[12]  discusses  NASTRAN  calculations  and  accompanying  tests  on  a  series  of 
integrally  stiffened  Wagner  beams.  The  results  of  the  numerical  calculations 
(stresses  at  the  crests  of  the  waves)  were  utilized  in  order  to  obtain 
equivalent  notch  stresses.  The  latter  formed  the  basis  of  a  fatigue  analysis 
which  was  correlated  with  the  experiments.  Correlation  factors  were  intro¬ 
duced  as  necessary.  The  main  emphasis  in  the  series  was  on  a  range  of  ulti¬ 
mate-shear/buckling  -  shear  ratios  of  2-3  (which  is  lower  than  the  range  to 
be  represented  in  this  report) . 

Reference  [13]  presents  a  series  of  fatigue  tests  on  stiffened  com 
posite  shear  panels.  Both  web  cracking  due  to  severe  bending  and  stiffened- 
web  separations  were  noted.  A  supporting  analysis  using  NASTRAN  was  also 
performed.  Interaction  of  the  buckled  waves  with  the  surrounding  structure 
were  not  taken  into  account  in  the  analysis.  Additional  tests  on  buckled 
composite  shear  panels  are  reported  in  [14,  15,  16,  19]. 


The  initiation  of  fatigue  cracks  at  the  points  of  the  most  severe 
bending  in  buckled  shear  panels  was  noted  in  Reference  [17]  .  The  tests 
were  conducted  on  aluminum  panels  and  the  results  agreed  reasonably  well 
with  estimates  based  on  SN  diagrams. 

Analytical  and  numerical  calculations  of  the  postbuckling  behavior 
of  stiffened  shear  panels  can  be  useful  for  the  prediction  of  the  crack 
nucleation  stage  of  the  fatigue  life.  Even  so,  they  are  rather  limited 
in  scope,  since  the  local  modification  of  the  stress  fields  due  to  inela¬ 
stic  strains  or  the  severe  localized  stress  concentrations  are  extremely 
difficult  to  obtain. 

The  analysis  can  predict  the  overall  behavior  including  the  number, 
amplitude  and  inclination  of  diagonal  waves  and  also  the  effects  of  the 
surrounding  structure  on  the  deformation  parameters.  Energy  methods 
appear  to  be  most  useful  in  this  respect  [18,  19] . 

The  numerical  calculations  can  use  one  of  several  computer  codes 
which  have  nonlinear  capability  for  plate  and  shell  elements.  The  prog¬ 
rams  which  were  most  extensively  used  were  NASTRAN  [20]  and  STAGS  [21] , 
Some  of  the  calculations  were  used  for  comparisons  with  experiments-  See 
for  example  [8,  13,  17,22,  23,  24].  It  appears  from  the  results  of  the 
calculations  that  overall  behavior  can  be  predicted  rather  accurately  but 
the  results  of  detailed  stress  distributions  are  less  conclusive. 

The  effects  of  the  surrounding  structure  (stiffeners,  stress  envi¬ 
ronments)  on  the  durability  of  repeatedly  buckled  panels  are  very  impor¬ 
tant.  Not  only  do  the  stiffeners  influence  the  overall  behavior,  but 


they  also  contribute  very  strongly  to  the  severe  local  effects  near  the 
corners  of  the  panel  -  where  most  fatigue  cracks  initiate.  In  fact,  many  of 
of  the  fatigue  failures  in  the  tests  reported  above  and  in  other  tests,  star¬ 
ted  at  the  stiffener  (or  flange)  panel  boundaries  near  the  corners. of  the 
tension  diagonal.  Thus  in  order  to  obtain  realistic  and  meaningful  information 
on  the  durability  of  shear  panels,  they  have  to  be  tested  in  a  realistic  envi¬ 
ronment,  with  realistic  stiffeners  and  flanges,  and  with  the  correct  transmi¬ 
ssion  of  the  shearing  loads  from  the  neighboring  structure. 

Based  on  the  above  considerations  and  on  preliminary  testing  Ill,  a 
three  point  loading  Wagner  beam  configuration  was  decided  upon  for  this  test 
series.  The  need  for  mutual  interaction  between  the  shear  panel  and  the  fle¬ 
xible  stiffeners  and  flanges  ruled  out  the  possibility  of  a  "picture  frame" 
testing  (the  effects  of  additional  moments  which  are  inherent  to  the  three- 
point  loading  configuration,  although  of  no  major  importance,  were  considered 
to  contribute  to  the  realistic  environment). 

The  stiffeners  and  flanges  of  the  Wagner  beams  were  bonded  to  the  shear 
panel  according  to  aeronautical  specifications.  The  fact  that  there  was  no 
bonding  fatigue  failure  during  the  entire  test  series  is  indicative  of  the 
integrity  of  the  bonding.  It  provided  a  realistic  support  condition  and,  on 
the  other  hand,  it  was  not  a  "weak  link"  in  the  test  system.  The  success  of 
the  bonding  in  this  series  of  rather  severe  postbuckling  cycling  tests  is  an 
important  result  although  it  was  not  an  initial  requirement  of  this  research 

It  is  recognized  that  the  durability  of  repeatedly  buckled  panels 
in  a  realistic  environment  is  both  a  multiparameter  problem  and  one  for  which 


the  statistical  nature  of  fatigue  behavior  should  apply.  Therefore,  any 
single  test  series  must  limit  its  scope  and  vary  only  a  few  of  the  parame¬ 
ters.  Within  the  Wagner  beam  configuration  which  was  chosen  to  be  tested, 
the  following  effects  were  studied  and  corresponding  parameters  varied: 

a.  Effect  of  the  surrounding  structure  -  by  varying  the  size  of  the  stiffness 
and  flanges  (three  sizes  were  used) . 

b.  Effects  of  the  magnitude  of  the  initial  buckling  load  -  by  varying  the  web 
thickness  (three  web  thicknesses  were  used). 


c.  Effects  of  the  shear  loading  (compared  with  the  initial  buckling  shear  a 
the  ultimate  shear  V^^)  -  by  varying  the  magnitudes  of  the  repeated  loads  Vc^c 
from  test  to  test.  The  range  of  the  applied  repeated  loans  was  on  the  high  side 
several  times  the  initial  buckling  loads  and  up  to  1/2  ^  .  The  main  conside¬ 

rations  for  this  choice  were  the  need  to  investigate  the  effects  of  severe  buck¬ 
ling  on  one  hand  and  to  maintain  the  loads  within  range  of  possible  fatigue 
spectra  of  existing  modern  day  airplanes  on  the  other  hand. 


A  constant  load  amplitude  of  0  -  V  was  maintained  within  each  test. 

1  eye 

Spectrum  testing  was  ruled  out  for  this  series  in  order  to  limit  the  number 
of  test  variables.  The  possibility  of  doing  some  spectrum  testing  in  a  fu¬ 
ture  series  is  envigased.  Such  results  could  be  checked  against  test-ori¬ 
ented  cumulative  damage  predictions  based  on  the  present  series. 

Material  properties  were  not  varied  in  this  test  series  and  all  the  tests 


were  conducted  on  2024T3  aluminum  specimens.  Extension  to  composite  beams  of 
similar  configurations  is  planned. 


Further  details  of  the  test  specimens,  setup  and  experimental  pro¬ 
cedures  are  given  in  the  subsequent  sections. 

Some  of  the  important  test  parameters  are  grouped  into  nondimen- 
sional  forms  and  a  correlation  study  was  initiated  in  order  to  arrive  at 
a  semi-empirical  "test  formula".  This  formula  should  be  checked  for  its 
predictive  capability  against  the  various  tests  performed  as  well  as 
against  future  tests.  Apart  from  its  obvious  practical  potential  as  a 
design  tool,  it  can  contribute  towards  the  understanding  of  the  relative 
importance  of  the  various  parameters  in  determining  the  fatigue  behavior. 

This  phase  of  the  research  is  still  under  investigation  and  will  be  con¬ 
tinued  throughout  the  entire  test  series,  including  future  tests.  Pro¬ 
gress  in  this  important  aspect  of  the  research  is  reported  in  a  subse¬ 
quent  section. 

The  STAGS  computer  program  was  chosen  for  numerical  support  to  the 
research  and  calculations  of  the  strong  postbuckling  behavior  of  the  pa¬ 
nels  and  surrounding  structure  were  conducted.  After  overcoming  initial 
difficulties,  good  correlation  between  tests  and  calculations  was  obtained 
by  taking  the  initial  imperfections  as  two  trigonometric  double  sine 
terms  with  amplitude  equal  to  the  web  thickness.  Apart  from  overall 
support  to  the  testing  program,  the  calculations  can  be  used  for  specific 
purposes  such  as  the  numerical  determination  of  the  load  level  at  which 
inelastic  stresses  first  occur  in  the  susceptible  region  near  the  corner. 

This  load  may  become  an  important  design  parameter  in  view  of  the  well  known 
influence  of  local  inelastic  stresses  on  low  cycle  fatigue  behavior.  Some  of 
the  numerical  results  are  presented  in  this  report. 


2.  TEST  SETUP  §  PROCEDURES. 


The  test  setup  is  shown  in  Fig.  1.  The  load  experienced  by  the 
beam  specimen  in  the  loading  frame  of  this  setup  is  of  a  three  point 
loading  type;  the  load  is  introduced  through  a  pin  into  the  center  of 
the  beam  by  a  hydraulic  jack  via  a  load  cell  connected  to  two  heavy 
parallel  vertical  bars,  into  which  the  pin  is  inserted.  The  load  is 
controlled  through  an  MTS  system  that  permits  precise  control  of 
applied  load  and  frequency  loading  ranges  from  zero  to  maximum  load 
(in  a  const,  amp.  mode).  The  edges  of  the  beam  are  loaded  by  pins  in¬ 
serted  into  the  specimen  and  into  two  heavy  parallel  vertical  bars  at 
each  end,  which  are  simply  connected  by  pins  to  the  loading  frame.  The 
reaction  loads  experienced  by  these  pins  are  introduced  via  the  uprights 
into  the  shear  web. 

Since  allowances  for  lateral  deflections  of  the  specimens  may  cause 
premature  failure  due  to  lateral  instability  of  the  flanges  of  some  of 
the  beam  configurations,  transverse  deflections  of  the  specimens  are 
prevented  by  means  of  heavy  parallel  cross  members  attached  to  the  loa¬ 
ding  frame  along  the  top  and  bottom  flanges  of  the  beam.  These  cross 
beams  permit  the  restraining  of  the  transverse  deflections,  only  at 
specified  points  along  the  specimen. 

The  response  of  the  specimens  is  measured  by  strain  gages  t  nded  to 
the  shear  web  surfaces,  flanges  and  uprights  of  the  beam  specimen  and  re¬ 
corded  by  a  multichannel  data  logger.  Recording  of  these  gages  is  li¬ 
mited  to  quite  a  low  number  of  cycles  relative  to  the  overall  life  of  the 


specimens,  because  the  strain  gages  cannot  sustain  the  very  high  cyc¬ 
ling  strains  experienced  by  the  beams  and  have  a  very  low  endurance 
life  under  such  circumstances. 

Also,  the  shadow-moire  technique  is  employed  for  overall  observa¬ 
tion,  of  the  progressive  behavior  of  the  buckled  shear  web,  as  well  as  for 
comparison  of  the  deflected  shape  of  the  beam  with  predictions  by  the 
STAGS  [21],  computer  code.  The  grids  employed  in  this  technique  are 
removed  when  cracks  are  assumed  to  initiate,  to  permit  a  visual  obser¬ 
vation  of  these  cracks  and  to  avoid  their  being  damaged  at  the  failure 
of  the  beam. 

Each  test  was  conducted  and  divided  into  the  following  phases: 

1.  A  slow  gradual  loading  of  the  beam  to  a  load  slightly  above  either 
the  predicted  or  observed  shear  buckling  of  the  specimen,  to  deter¬ 
mine  the  actual  buckling  load  of  the  webs  of  the  beam.  Then  release 
of  the  load. 

2.  A  slow  gradual  reloading  of  the  beam  to  the  load  of  (1)  to  identify 
any  effects  of  the  load  in  (1)  on  the  critical  load  observed  in  (2) . 
Then  either  release  of  the  load,  or  further  increase  of  the  load. 

3.  In  this  phase  two  different  procedures  of  loading  were  practiced;  In 
the  first  one,  gradual  loadings  of  the  specimen(WB-t)to  determine  the 
"permanent  buckling"  load,  was  performed.  The  load  was  increased 

by  specified  incraments  relative  to  the  preceding  load  and  then  re¬ 
leased.  This  load  was  gradually  increased  in  this  fashion  up  to  the 
maximum  load  for  repeated  loading,  the  value  of  which  was  dictated 
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by  either  the  strain  value  at  the  lower  comer  of  the  tension  dia¬ 
gonal,  or  the  planned  ratio  of  the  repeated  working  load  to  the 
predicted  ultimate  load.  In  some  cases  the  magnitude  of  the  working 
load  was  dictated  by  both  the  measured  strain  at  the  corner  and  the 
planned  loading  ratio.  It  should  be  noted  that  the  strain  value  at 
the  lower  corner  is  usually  well  in  the  plastic  region  and  sometimes  very 
high  strains  may  be  observed  prior  to  reaching  the  planned  loading  ra¬ 
tio.  Therefore,  further  loading  of  the  specimen  can  result  in  very 
high  strains,  which  will  trigger  failure  of  the  beam  after  a  very  short 
cycling  life.  Consequently,  repeated  loading  is  executed  at  a  smaller 
loading  ratio  than  the  planned  one. 

In  the  second  procedure (Beams  WB-2  to  WB-5  and  WB-7  to  WB-11),  the  load 
is  continuously  increases  to  the  maximum  load  for  repeated  loading. 

The  magnitude  of  this  load  is  determined  as  in  the  first  procedure. 

4.  Continuous  repeated  loading  from  V=0  to  with  interruptions  every 

few  tousand  cycles  for  complete  strain  gage  and  moire  records.  The 
beams  are  loaded  at  very  low  frequency  (1  to  2.2  Hz). 

5.  Once  a  crack  initiates  and  completely  penetrates  through  the  thick¬ 
ness  of  the  web,  its  further  growth  is  continuously  measured,  moni¬ 
tored  and  recorded  by  a  video  camera  and  tape  till  large  scale  fai¬ 
lure  of  the  web  takes  place.  It  should  be  noted  that  in  order  to 
avoid  complete  failure  of  the  whole  beam,  due  to  failure  of  the  web, 
a  stop  is  introduced  in  the  center  of  the  beam,  thus  avoiding  any 
significant  damage  to  the  beam  at  collapse. 


ft*** 
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5.  TEST  SPECIMENS 


10  Shear  beams  were  tested.  Their  dimensions  are  presented  in  Fig. 2 
and  Table  1.  It  can  be  seen  that  the  test  series  comprises  of  beams  with 
three  web  thicknesses:  0.5  mm,  0.8  mm  and  1.0  mm  and  stiffened  by  three 
types  of  L  sections:  "light"  -  112*12x1 .0  mm,  "moderate"  -  L15xl5*1.5  mm 
and  "heavy"  -  L20*20*2  mm.  This  allows  for  evaluation  of  the  effects  of 
quite  a  wide  range  of  stiffening  and  flexibility  of  the  framing,  as  reflec¬ 
ted  by  the  parameter  (I^./b3t)  and  of  the  corresponding  relative  magnitude 
of  the  buckling,  yield  and  ultimate  loads  of  the  beams  (see  Table  2)  on  the 
durability  of  the  beam  at  the  specified  cyclic  loading  condition  (see  Table  2). 

The  critical  shear  and  ultimate  loads  corresponding  to  the  beams 
given  in  Table  2,  were  predicted  by  Eqs.  (32)  &  (33a)  of  [3]  respectively, 
as  well  as  by  the  STAGS  computer  code  [21].  Since  this  program  permits  a 
comprehensive  detailed  analysis  of  the  beam,  it  provides  information  about 
the  "critical"  points  along  and  adjacent  to  the  tension  diagonal  of  the  web 
and  initiation  of  yielding  there.  This  information  is  compared  and  corre¬ 
lated  later  with  the  recorded  data  by  the  strain  gages  adjacent  to  these 
"critical"  points,  to  obtain  the  actual  load  at  initial  yielding  of  the  beam. 

As  can  be  seen  in  Fig.  2 ,  the  beams  consist  of  five  shear  webs  sym¬ 
metrically  framed  by  identical  L  section  flanges  and  uprights.  Two  of  the 
webs  serve  as  the  test  section  fields.  Through  the  other  three,  the  central 
one  and  the  edge  ones,  the  three  loads  are  introduced  into  the  beam.  These 
webs  are  heavily  stiffened  in  a  manner  which  assures  almost  uniform  shear 
load  difusion  from  the  three  loading  pins  into  the  test  panels. 
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In  the  early  stages  of  the  test  program  (beams  WB-2  to  WB-7) ,  the 
two  test  sections  were  tested  simultaneously.  This  was  aimed  at  actually 
doubling  the  number  of  test  data  points,  as  well  as  confirming  the  almost 
symmetrical  loading  behavior  of  the  test  rig.  However,  the  actual  ultimate 
load  was  required  for  each  specimen  to  determine  its  actual  cycling  loading 
ratio  and  therefore  it  was  decided  to  use  each  beam  both  for  cycling  and 
static  testing  (of  course  after  symmetry  and  unbiased  loading  of  the  rig 
was  evident  from  the  first  tests) .  This  was  carried  out  by  heavily  stif¬ 
fening  one  of  the  test  sections,  so  that  it  would  be  damaged  as  little  as 
possible  during  the  cycling  loading  phase  up  to  the  total  collapse  of  the 
other  test  section.  Then,  this  section  was  utilized  for  the  ultimate 
static  test  in  the  inverted  manner  described  earlier.  3  pairs  of  strain  gage  ro 
settes  were  bonded  face  to  face  on  both  sides  of  the  shear  web.  They  were 
located  at  the  two  "critical"  zones  at  each  end  of  the  tension  diagonal 
as  close  as  possible  to  the  frame,  and  at  the  center  of  the  panei(Fig.  2).  This 
latter  one  was  used  for  reference  measurements,  as  well  as  for  determining 
the  buckling  load  and  overall  average  behavior  of  the  beam. 


4.  NUMERICAL  STUDIES 


The  STAGS  computer  program  [21]  was  chosen  for  numerical  support  of 
the  research. 

Calculations  of  the  buckling  load  and  postbuckling  behavior  of  the 
panels  and  surrounding  structure  were  conducted  for  beams  WB-3  through  WB-8 
and  WB-11.  The  analytical  model  used  to  represent  the  tested  field  of  beams 
WB-3,  WB-4  and  WB-7  was  a  160x160x0.5  mm  plate  attached  to  a  stiffening  frame 
comprising  of  a  thicker  plate  with  a  width  of  20  mm  and  thickness  of  4.5  mm. 
(The  need  to  replace  the  L  stiffeners  by  a  thicker  plate  stems  from  the  ina¬ 
bility  of  the  STAGS  version  used  in  the  calculation  to  account  for  the  in¬ 
plane  rigidity  of  the  stiffeners) . 

For  beams  WB-5  and  WB-11,  the  dimensions  of  the  plate  were  160x160x0.8  mm 
and  the  thickness  of  the  frame  was  4.8  mm,  and  for  beams  WB-6  and  WB-8,  the 
dimensions  of  the  plate  were  160x160x1.0  mm  and  the  thickness  of  the  frame 
was  5.0  mm. 

The  plate  was  attached  to  elastic  beams  which  simulated  the  properties 
of  the  vertical  stiffeners  (r  eights)  and  flanges  of  the  tested  beams  (see 
Fig .  3  and  Table  3) . 

The  vertical  stiffeners  were  loaded  uniformly  by  distributed  shear. 

In  addition,  forces  normal  to  the  cross  section  of  the  flanges,  were  introduced 
to  simulate  the  bending  moment  acting  on  the  beam. 

The  plate  was  discretized  into  a  29x29  non  uniform  grid  with  denser  meshes 


at  the  corners. 


The  numerical  analysis  with  the  STAGS  was  performed  in  the  two  follo¬ 
wing  modes : 

a.  A  bifurcation  analysis  to  predict  the  buckling  load  of  the  shear  webs 
of  the  beam.  The  initial  buckling  load  is  essential  to  the  analysis 
of  the  repeatedly  buckled  beams. 

b.  An  elastic  nonlinear  postbuckling  analysis  with  assumed  initial  im¬ 
perfections  to  support  and  supplement  the  test  results.  An  important  ob¬ 
jective  of  this  phase  of  the  numerical  analysis  was  also  the  definition 
of  the  shear  load  where  yielding  first  occurs  at  the  critical  lo¬ 
cations  near  the  tension  diagonal. 

The  buckling  load  yielded  by  the  numerical  studies  with  STAGS  was  com¬ 
pared  with  the  calculated  buckling  load  obtained  with  Eq  32  of  Ref.  3  and 
with  the  actual  buckling  load  observed  at  the  test.  These  critical  loads 
are  given  in  Table  2. 

In  the  early  studies  with  STAGS,  the  imperfection  w^  of  the  web  was 
taken  as  a  one  term  half  sine  wave  in  the  x  and  y  directions. 

wQ(x,y)  =  aQsin  (-^  j  sin(^) 

The  choice  of  imperfection  mode  was  for  the  purpose  of  observing  the  deve¬ 
lopment  and  emergence  of  changing  buckling  and  post-buckling  deformation 
patterns.  Application  of  this  imperfection  mode  has  indicated  that  severe 
convergence  problems  were  encountered  unless  an  imperfection  amplitude  a^ 
of  a  magnitude  equal  to  the  thickness  of  the  web  was  introduced  into  the 
calculations.  Also,  increasing  experience  with  the  above  one  term  imper¬ 
fection,  w^Cx.y),  proved  it  to  be  insufficient  for  studying  the  post  buckling 


behavior  of  the  beam.  Therefore,  the  following  two  term  expression 


wQ(x,y)  =  a( 
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was  adopted  for  further  studies  with  the  STAGS  program. 


In  the  numerical  studies,  the  shear  load  V  was  increased  in  incre¬ 
ments  up  to  a  maximum  shear  which  was  well  over  the  working  cyclic  loa¬ 
ding.  The  maximum  load  for  WB-3,  WB-4  and  WB-7  was  1400  kg  and  for 
WB-5  and  WB-11,  the  load  was  1500  kg.  No  convergence  problems  were  en¬ 
countered  in  this  range. 


The  computed  load-vs-def lections  curves  as  well  as  the  corresponding 
experimental  curves  are  shown  in  Fig.  4  and  Fig.  5  for  WB-3  and  WB-5  res¬ 
pectively.  The  somewhat  higher  deformations  obtained  in  the  numerical  cal 
culations  are  probably  due  to  the  more  severe  imperfection  used  in  the 
calculations  (compared  with  the  experiments)  and  slightly  more  flexible 
support  conditions.  Apart  from  this,  the  numerical  results  follow  the 
test  results  quite  well. 


In  Figs.  6  and  7  the  calculated  deflected  shapes  along  the  two 
diagonals  of  the  web  of  beam  WB-3  at  different  load  levels,  are  compared 
with  the  test  results.  A  good  correlation  between  the  predicted  shapes 

and  the  experimental  results,  is  observed.  The  scatter  at  some  points 

* 

is  due  to  difficulties  in  analysing  the  moire  records  in  particular  in 
the  small  deflection  range.  Similar  results  were  obtained  for  the  other 
tested  beams. 


The  deflected  pattern  of  beam  WB-7  obtained  from  the  STAGS  program 
for  V  =  600  kg  is  presented  in  Fig.  8  and  compared  with  the  corresponding 
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5.  TEST  RESULTS  -  DISCUSSION  §  ANALYSIS. 

The  test  results  are  presented  in  Table  2,  together  with  the  pre¬ 
dicted  results  using  the  numerical  studies  with  STAGS  [21]  and  those 
calculated  according  to  Ref.  [3  . 

a.  Buckling  loads. 

The  experimental  buckling  loads  (^cr)eXp»  were  obtained  from  load- 
strain  responses,  recorded  by  a  pair  of  strain  gages  bonded  face  to  face 
at  the  center  of  the  web  (see  Section  3).  The  critical  loads  were  deter¬ 
mined  by  the  method  described  in  Fig.  9  for  beam  WB-11.  It  is  seen  from 
thi  s  table  that  a  rather  good  agreement  exists  between  the  calculated  loads 
by  the  two  methods  of  analysis  and  that  correlation  between  the  theoreti¬ 
cal  predictions  and  the  experimental  buckling  load  depend  upon  the  relative 
rigidity  of  the  frames. 

Fair  agreement  is  observed  for  the  beams  with  the  "heaviest"  frames, 
WB-3,  WB-4  and  WB-7,  whereas  very  good  correlation  is  found  for  the  beams 
with  the  relatively  moderate  stiffening,  WB-5,  WB-6,  WB-8  and  WB-11,  and 
the  most  flexible  frames  WB-9  and  WB-2.  This  implies  that  the  actual 
rigidity  of  the  stiffeners  in  the  experiments  is  less  than  that  assumed 
in  the  calculations.  The  effect  is  more  pronounced  in  the  case  of  the 
heavier  stiffeners,  which  (theoretically)  should  have  approached  the 
case  of  fixed  supports. 


A  typical  moire  fringe  picture  of  a  buckled  web  showing  the  deflection 
pattern  at  the  buckling  load  is  given  in  Fig.  10  for  beam  WB-6.  The  shape 


of  the  buckling  mode  is  in  agreement  with  the  results  of  the  STAGS  cal¬ 
culations.  It  can  be  useful  for  further  improved  calculations. 

b.  Ultimate  failure  load. 

In  the  static  ultimate  tests  it  was  observed  that  collapse  of  the  beams 
was  associated  with  two  different  modes  of  failure:  Shearing  of  the  web 
along  a  diagonal  perpendicular  to  the  tension  field  in  beams  WB-3,  WB-4, 
WB-5  (see  Fig.  11),  or  failure  of  the  flanges  of  the  beam,  either  in 
yielding  or  buckling^for  the  rest  of  the  beams  (see  Fig.  12).  It  is  seen 
from  Table  2,  that  neither  of  these  modes  of  failure  had  a  definite  effect 
on  the  agreement  between  the  predicted  ultimate  loads  (eq.2  of  Ref. [3])  and 
the  experimental  ones.  However,  the  correlation  for  the  beams  which  failed 
in  the  "tearing  mode",  resulted  in  a  ratio  of  the  experimental  load  to 
the  calculated  load  which  was  equal  to,  or  greater  than  unity,  whereas 
the  ratios  for  the  other  beams  varied  between  values  smaller  than  unity 
and  above  unity. 

It  should  be  noted  that  the  flanges  of  beams  WB-6  and  WB-8  experienced 
local  bucklings  at  loads  which  were  considerably  below  the  predicted  ul¬ 
timate  loads  (see  Table  2).  After  this  behavior  was  observed  for  beam 
WB-6  (and  later  for  beam  WB-8),  loading  of  beam  WB-8  was  stopped  immedia¬ 
tely  at  incipient  buckling.  The  beam  was  unloaded,  its  compression  flange 
was  stiffened  and  the  beam  was  subsequently  reloaded  untill  collapse. 

It  is  apparent  from  Table  2  that  stiffening  of  the  flange  had  a  signifi¬ 
cant  effect  on  the  failure  load  of  the  beam  and  resulted  in  a  load  almost 


equal  to  the  predicted  one. 


c.  Repeated  bucklings. 

Table  2  presents  the  number  of  cycles  to  initiation  of  a  crack  at 
the  critical  comer  (as  actually  first  observed  during  the  test;  see 
later  discussion) ,  and  the  number  of  cycles  to  complete  failure  of 
the  shear  web  for  the  prescribed  cyclic  loading  conditions  represented 

^ult^  cal  .^cyclic  anc*  ^cyclic^cr  ‘ 

During  the  repeated  buckling  tests,  cracks  were  first  detected,  as 
expected^at  the  lower  corner  of  the  web  along  the  tension  diagonal,  ex¬ 
cept  for  beam  WB-11  where  the  crack  was  first  observed  at  the  upper 
corner  of  the  diagonal.  These  cracks  whre,  when  observed,  of  a  length 
of  about  ten  times  the  thickness  of  the  web  and  it  is  reasonable  to 
assume  that  they  were  formed  by  the  severe  bending  stresses  on  the  ten¬ 
sion  bending  side  at  the  corner. 

Since  these  cracks  formed  at  the  very  edge  of  the  web  almost  under 
the  stiffeners  and  they  first  developed  by  propagation  through  the  thick¬ 
ness  of  the  web  from  the  tension  to  the  compression  side,  their  detec¬ 
tion  was  almost  impossible  before  they  have  reached  a  certain  length  and 
depth.  Therefore,  the  value  of  should  be  regarded  with  caution.  Also 
it  was  observed  that  at  this  stage,  the  cracks  mainly  propagated  through 
the  thickness  of  the  web  and  were  almost  "steady"  lengthwise.  After  com¬ 
prehensive  penetration  through  the  thickness,  the  cracks  began  growing 
lengthwise  and  their  growth  was  measured  (see  Fig.  13).  The  cracks  star¬ 
ted  either  horizontally  and  immediately  followed  by  a  vertical  branch 


or  vice  versa. 


As  can  be  seen  from  Table  2  the  rate  of  growth  of  the  crack  and  the 


marginal  cyclic  life  of  the  beam  to  reach  depended  upon  the  working 

load  ratio,  V  , •  /V  yielding  low  N..-N.  values  for  relatively 

*  cyclic  ult’  f  i 

high  load  ratios.  It  was  also  observed  that  in  the  case  of  low  load  ra¬ 
tios  the  cracks  propagated  lengthwise  in  a  non  continuous  fashion  up  to 
complete  failure,  whereas  for  the  higher  ratios,  the  behavior  was  quite 
smooth. 


d.  Correlation  studies. 


The  aim  of  the  present  work  is  to  obtain  design  charts  for  predictions 
of  durability  of  post  buckled  shear  webs  and  design  recommendations. 
Therefore,  an  attempt  was  made  to  establish  a  design  curve  from  the  pre¬ 
sent  test  results  and  the  corresponding  analytical  studies.  It  has  been 
tentatively  assumed  that  the  durability  of  the  shear  web  is  mainly  go¬ 


verned  by  the  following  nondimensional  parameters:  vCyCiic/v 
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It  should  be  noted  that  based  on  the  discussion  regarding  the  accu¬ 
racy  in  in  the  previous  paragraph,  it  seems  to  be  more  meaningful  to 
predict  rather  than  N^, 

A  least-square  method  was  used  to  obtain  the  values  of  A,  a,  g,  y, 
which  resulted  in:  A  =  0.81,  a  =  -0.81,  8  =  4.91,  y  =  0.86. 


It  is  seen  from  Table  2  that  Nf  for  beam  WB-3  deviates  significantly 


from  the  value  that  might  be  expected  based  on  the  values  experienced 
for  the  identical  beams  WB-4  and  WB-7.  Note  that  beams  WB-4  and  WB-7 
yielded  =  22300  at  a  working  load  of  900  kg  and  =  168410  at  a 
working  load  of  600  kg  respectively.  The  corresponding  values  for  WB-3 
were  =  23000  at  700  kg.  This  last  result  for  appears  to  be  un¬ 
reasonable,  since  a  much  higher  value  was  expected.  The  possibility  of 
a  defective  beam  should  be  considered  and  a  retest  is  scheduled.  This 
deviation  also  affected  significantly  the  least  square  results,  when  it 
was  included  in  the  studies.  Therefore  it  was  excluded  from  the  analysis. 
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6.  CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FUTURE  WORK. 


1.  At  load  levels  within  the  flight  envelope  (27%  to  46%  of  ultimate 
loads) ,  repeated  buckling  led  to  initiation  (N^)  and  propagation  of 
fatigue  cracks,  with  subsequent  failure  (N^.) .  Typical  values  of 
numbers  of  load  cycles  N^  and  N^  ranged  from  about  20^000  cycles  to 
200,000  cycles.  These  values  demonstrate  the  seriousness  of  the  prob¬ 
lem  to  the  current  design  of  flight  vehicles. 

2.  Cracks  started  at  the  most  highly  stressed  regions  (as  verified  by  the 
numerical  analysis  and  by  observed  curvatures  of  the  wave  patterns) 

at  the  corners  near  the  tension  diagonals,  along  the  attachment  lines 
to  the  stiffeners  and  flanges. 

3.  A  design  formula  was  developed,  which  can  predict  logN^  in  this  test  se¬ 
ries  to  within  a  scatter  of  19%. 

4.  From  the  formula,  it  appears  that  the  most  significant  parameter  affec¬ 
ting  N,.  was  the  ratio  of  V  l4/V  .  The  effects  of  the  other  para- 

&  f  ult  cyclic  v 

meters  were  close  to  linear. 

5.  An  increase  in  the  relative  rigidity  of  the  frames  reduced  the  life  N^ 
of  the  structure.  Thus,  it  seems  that  over-rigidity  of  the  supporting 
frames  near  the  corners  may  lead  to  premature  failure. 

6.  Good  agreement  was  obtained  between  the  test  results  and  the  supporting 
numerical  analysis.  This  can  be  seen  from  comparisons  of  the  deflection 
patterns  of  the  web, (in  Figs.  6,  7,  and  8). 


Based  on  the  current  results  and  conclusions,  investigations  are  being 

persued  in  the  following  directions: 

a.  Design  formulas  are  now  being  developed  for  predicting  crack  ini¬ 
tiation  (N^) ,  based  on  an  experimental  "standardized"  minimum  crack 
size. 

b.  A  similar  investigation  will  be  conducted  on  the  critical  crack  size 
and  on  the  rate  of  crack  propagation. 

c.  An  investigation  will  be  conducted  on  the  influence  of  additional 
parameters  and  their  predictive  capability.  An  important  possibi¬ 
lity  in  this  respect  which  is  now  under  study  is  V  ,  which  is  the 
shearing  load  at  which  local  yielding  first  starts  at  the  location 
of  crack  initiation. 

d.  Efforts  to  improve  the  numerical  and  analytical  support  capabilities 
are  now  being  made. 
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Table  No.  2  Predictions  and  Test  Results 
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Fig.  3  MODEL  FOR  NUMERICAL  STUDY  BY  STAGS  PROGRAM 
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Fig.  11  TYPICAL  "TEARING  FAILURE"  MODE 


Fig.  12  TYPICAL  FAILURE  MODE  IN  BENDING  AND  BUCKLING  OF  THE  FLANGE 


